Accurate determination of elastic parameters for multi-component membranes 
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Heterogeneities in the cell membrane due to coexisting lipid phases have been conjectured to 
play a major functional role in cell signaling and membrane trafficking. Thereby the material 
properties of multiphase systems, such as the line tension and the bending moduli, are crucially 
involved in the kinetics and the asymptotic behavior of phase separation. In this Letter we present 
a combined analytical and experimental approach to determine the properties of phase-separated 
vesicle systems. First we develop an analytical model for the vesicle shape of weakly budded biphasic 
vesicles. Subsequently experimental data on vesicle shape and membrane fluctuations are taken and 
compared to the model. The combined approach allows for a reproducible and reliable determination 
of the physical parameters of complex vesicle systems. The parameters obtained set limits for the 
size and stability of nanodomains in the plasma membrane of living cells. 

PACS numbers: 87.15.Ya, 87.16.Dg, 87.17.Aa, 02.40.Hw 
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The recent interest in coexisting phases in lipid bilay- 
ers originates in the supposed existence of lipid hetero- 
geneities in the plasma membrane of cells. A significant 
role in cell signaling and traffic is attributed to small lipid 
domains called "rafts" [l], ■ While their existence in liv- 
ing cells remains the subject of lively debate, micrometer- 
sized domains are readily reconstituted in giant unilamel- 
lar vesicles (GUVs) made from binary or ternary lipid 
mixtures [3( . Extensive studies of these and similar model 
systems have brought to light a rich variety of phases, 
phase transitions and coexistence regimes [J|. In con- 
trast to these model systems, no large (micrometer sized) 
membrane domains have been observed in vivo. If in- 
deed phase separation occurs in vivo, additional processes 
which can arrest it prematurely must be considered. It 
has been suggested that nanodomains might be stabi- 
lized by entropy @ or that, alternatively, active cellular 
processes are necessary to control the domain size [f|. 
A third explanation is that curvature-mediated interac- 
tions might conspire to create an effective repulsion be- 
tween domains, impeding and ultimately halting their 
fusion as the phase separation progresses. Each of these 
three processes depends critically on membrane parame- 
ters such as line tension Q, curvature moduli and even 
the elusive Gaussian rigidities §]. Although some stud- 
ies report values [9( or upper bounds 0, H3 for these 
membrane parameters, a systematic method to deter- 
mine them from experiments that does not require ex- 
tensive numerical simulation and fitting is lacking. We 
present here a straightforward fully analytical method 
that allows for a precise, simultaneous determination of 
the line tension, the bending rigidity and the difference 
in Gaussian moduli from biphasic GUVs. Both the liquid 
ordered L a and the liquid disordered phase are quan- 
titatively characterized with high accuracy. Our method 
relies on an analytical expression for the shape of a mod- 




FIG. 1: Fluorescence raw data (red: L domain, green: L^ 
domain) with superimposed contour (light blue) . Insets: prin- 
ciple of contour fitting; a: intensity profile normal to the vesi- 
cle contour (taken along the dashed line in the main image); 
b: first derivative of the profile with linear fit around the vesi- 
cle edge (white line). The red point marks the vesicle edge. 

erately budded vesicle. A one-parameter fit to experi- 
mental shapes permits unambiguous determination of the 
line tension and the difference in Gaussian moduli. Our 
results provide important clues as to the origin and mag- 
nitude of long-ranged membrane- mediated interactions, 
which have been proposed recently as an explanation for 
the trapped coarsening 12, 13| and the very regular do- 
main structure of a meta-stable state 1J] found in exper- 
iments. Furthermore, our results show that nanometer- 
sized phase separated domains will be stable in life cells. 

Model The free energy associated with the bending 
of a thin membrane is described by the Canham-Hclfrich 
free energy |L5[. We ignore any spontaneous curvature 
of the membrane because the experimental system has 
ample time to relax any asymmetries between the leaflets. 
For a two-component vesicle with line tension r between 
the components, the free energy then reads: 
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where the Ki and the bending and Gaussian mod- 

uli of the two phases, respectively, the <Ji are their surface 
tensions, and p is the internal pressure. In equilibrated 
shapes such as our experimental vesicles, the force of the 
internal Laplace pressure is compensated by the surface 
tensions; consequently, both contributions drop out of 
the shape equations [161 ]. For each phase, we integrate 
the mean (H) and Gaussian (K) curvature over the mem- 
brane patch Si occupied by that phase; the line tension 
contributes at the boundary dS of the two phases. Using 
the Gauss-Bonnet theorem, we find that the Gaussian 
curvature term yields a constant bulk contribution plus 
a boundary term [tH |. 

The axisymmetric shapes of interest (Fig. [TJ are fully 
described by the contact angle ip as a function of the 
arc length s along the surface contour. The coordi- 
nates (r(s), z(s)) are fixed by the geometrical conditions 
f = cos^(s) and z = — sinip(s), where dots denote 
derivatives with respect to the arclength. Variational cal- 
culus gives the basic shape equation [161 ]: 
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This equation holds for each of the phases separately. 
The radial coordinate r(s) and tangent angle ip(s) must 
of course be continuous at the domain boundary. Addi- 
tionally, the variational derivation of equation @ gives 
two more boundary conditions [l8j |: 

Iim(/C2^(e) - MV>(-£)) = {Ak + Ak g )^^, (3) 

\im(K 2 ip(e) — Kix/j(—e)) 

elO 

(oa , a ^cosi/'osinV'o . sinV>o , . 
= -(2Ak + Ak g ) ^ — T ' ( 4 ) 
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with R n and ipo the vesicle radius and tangent angle at 
the domain boundary, Ak = Ki — k 2 , Ak g = Kq 1 ' — Kq\ 
and the domain boundary located at s = 0. 

The sphere is a solution of the shape equation we 
can therefore use it as an ansatz for the vesicle shape 
far from the domain boundary. We split the vesicle into 
three parts: a neck domain around the domain bound- 
ary, where the boundary terms dominate the shape, and 
two bulk domains, where the solution asymptotically ap- 
proaches the sphere. Perturbation analysis, performed 
by expanding Eq. ^ around the spherical shape, then 
gives for the bulk domains: 



, ('0 
/(*) r \ s + s o 

Vw s ) = — 



c i R l log I — 



(5) 



Here Ri is the radius of curvature of the underlying 
sphere and Sq the distance (set by the area constraint on 



the vesicle) from the point r = to the domain boundary. 
As was shown by Lipowski Q , the governing length scale 
is the invagination length, defined as the ratio £j = Ki/r 
of the bending modulus and the line tension. Our three- 
domain approach applies when this invagination length 
is small compared to the size of the vesicle. At s = the 
line tension, rather than the bending modulus, becomes 
the dominant term in the energy. Self-consistency of the 
solution requires that the deviation from the sphere solu- 
tion at that point be small, i.e. given by the dimension- 
less quantity Ri. This fixes the integration constant q. 

Near the domain boundary, ip must have a local ex- 
tremum in each of the phases and we can therefore ex- 
pand it as 

(6) 



These neck solutions must match at the domain bound- 
ary and also satisfy conditions [3J . Furthermore, they 
also need to match the bulk solutions to ensure continuity 
of "0 and its derivative "0. In total this yields seven equa- 
tions for the eight unknowns {iPq\ ipo\ rfoK s i}- The 
necessary eighth equation is provided by the condition of 
continuity of r(s) at the domain boundary. 

Put together the neck and bulk components of ip give 
a vesicle solution for specified values of the material pa- 
rameters {Ki, Ak g , t}. This solution compares extremely 
well to numerically determined shapes (obtained using 
the Surface Evolver package [lS], Fig. [3]). Moreover, for 
the symmetric case of domains with identical values of k, 
we can compare to earlier modeling in Ref. Q . The vesi- 
cle can then be described by a single dimensionless pa- 
rameter A = £/i?o, where 47ri? 2 equals the vesicle area. 
The 'budding transition' (where the broad neck desta- 
bilizes in favor of a small neck) is numerically found in 
Ref. 0] to occur at A = 4.5 for equally sized domains; 
our model gives a value of A = 4.63. 

Experiment Giant unilamellar vesicles (GUVs) were 
produced by electroformation from a mixture of 30 % 
DOPC, 50 % brain sphingomyelin, and 20 % choles- 
terol at 55 °C. Subsequently lowering the temperature 
to 20 °C resulted in the spontaneous formation of liq- 
uid ordered L Q and liquid disordered domains on the 
vesicles. The phase was stained by a small amount 
of Rhodamine-DOPE (0.2 %). In order to stain the L Q 
phase a small amount (0.2 %) of the ganglioside GM1 
was added, and subsequently choleratoxin labeled with 
Alexa 647 was bound to the GM1 [20j] . For imaging we 
chose a wide-field epi-fluorescence setup because short 
illumination times (1-5 ms) prevent shape fluctuations 
with short correlation times from being washed out. The 
raw data of a typical vesicle is shown in Fig. [TJ The 
lateral resolution of the equatorial optical sections was 
limited by diffraction and pixelation effects. In the nor- 
mal direction, however, a high (sub-pixel) accuracy was 
obtained. The upper inset in Fig.[T]shows a typical inten- 
sity profile along a line perpendicular to the contour. We 
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determine numerically the profile's first derivative (lower 
inset in Fig. [1]) and fit the central part around the max- 
imum intensity with a straight line. The intercept with 
the x-axis gives the position of the vesicle edge. The posi- 
tional accuracy achieved is typically 20 run. The contours 
obtained were subsequently smoothed by a polynomial 
and all contours from the same vesicle (typically around 
I 000) were averaged to give the final result for the mean 
contour. 
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FIG. 2: Fluctuation spectra of the ordered (red circles) and 
disordered (green circles) domains. The corresponding best 
fits of Eq. are shown in blue and black respectively. Inset: 
Typical real-space fluctuations along the vesicle perimeter. 

Spectra of the shape fluctuations were obtained from 
those parts of the contours that were nearly circular, 
i.e. far away from the neck domain. Fluctuations were 
determined for each single contour as the difference be- 
tween the local radius r and the ensemble averaged ra- 
dius R of a circle fitted to patches around the vesi- 
cles' poles: u(s) = r(s) — R where s is the arclength 
along the circle, see Fig. [2] The experimental fluctu- 
ation spectrum was obtained by Fourier transform as 
Uk = a J-a/2^ s r{s)e~ lk ' s , where a is the arclength of 
the contour patch, and k — n ■ — with n a non-zero in- 
teger. Taking into account the finite patch size 2l| and 
following the spectral analysis of a closed vesicle shell de- 
veloped by Pecreaux et al. [22j leads to a power spectrum 
for the vesicle fluctuation 
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Here q — ^-m with m a non-zero integer, L — 2ttR, 
and (|M g | 2 ) sph the spectrum of the entire vesicle derived 
in [22j], where the overline indicates temporal averag- 
ing during the illumination time. Eq. {7} was derived 
from the Canham-Helffich free energy for a flat mem- 
brane with periodicity L. However, as shown in [23], the 
spectrum of a sphere with radius R differs from that of 
the flat membrane only for the lowest wave numbers k. 
Therefore we can use Eq. ([7]) to fit our fluctuation spec- 
tra if we omit the two lowest modes. Examples of such 
fits are shown in Fig. [2] 

Results Fits of the fluctuation spectra using Eq. ([7|) 
give the values of the bending moduli and surface ten- 
sions of the two phases. Using these values, we fit the 
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FIG. 3: Example for an experimentally obtained ip(s) plot 
(red: L phase, green: phase ) together with the best fit 
of the model (blue: L phase, black: phase). The dashed 
lines mark the transition points between the neck and bulk 
domains. Insets: Fit to a numerically obtained shape (using 
Surface Evolver). 

experimentally obtained vesicle shapes with the model 
described above. This leaves us with two parameters: 
the line tension t between the two phases and the dif- 
ference Ak g between their Gaussian moduli. Since the 
experimental data show that tp at the domain boundary 
follows a straight, continuous line we further assume that 
the derivative ip is continuous at the domain boundary 
(as suggested before 



1S(] ) . Imposing this additional 



condition fixes the value of Ak g for given r, leaving us 
with a single free parameter to fully describe the sys- 
tem [23|. By fitting the experimental data, we directly 
extract the line tension. An example fit is shown in Fig. 13; 
Values found for the bending moduli are 8± I • I0~ 19 J for 
the L Q domain and 1.9 ± 0.5 • I0~ 19 J for the Ld domain. 
For the line tension we found a value of 1.2 ± 0.3 pN, 
which is in the same range as that estimated by Baum- 
gart et al. ■ Finally, the difference in Gaussian moduli 
is about 3 ± I ■ 10~ 19 J, in accordance with the earlier es- 
tablished upper bound (k g < —0.83k) reported by Siegel 
and Kozlov hfl|. An overview of the results is given in 
Table H 

Discussion Ultimately, one worries about the mem- 
brane's elastic parameters because their precise magni- 
tude has important consequences for the morphology and 
dynamics of cells. The literature is replete with theoreti- 
cal speculations which depend strongly on, among others, 
the line tension. While the values we report apply to re- 
constituted vesicles, we can nonetheless use them in some 
of these models to explore possible implications for cellu- 
lar membranes. The majority of the investigated vesicles 
finally evolved into the fully phase separated state. This 
finding is in agreement with previous work by Frolov et 
al. [5j , which predicts, for line tensions larger than 0.4 pN, 
complete phase separation for systems in equilibrium. 
It should be noted that the line tension found is also 
smaller than the critical line tension leading to budding: 
recent results by Liu et al. 24| show that for endocytosis 
by means of membrane budding both high line tensions 
(> 10 pN) and large domains are necessary. Therefore 
nanodomains will be stable and will not bud. In cells, 
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however, additional mechanisms must be considered. To 
explain the absence of large domains in vivo, Turner 
et al. 0] make use of a continuous membrane recycling 
mechanism. For the membrane parameters we have de- 
termined such a mechanism predicts asymptotic domains 
of ~10 nm in diameter. Our results, in combination with 
active membrane recycling therefore support a minimal 
physical mechanism as a stabilizer for nanodomains in 
cells. A separate effect, purely based on the elastic prop- 
erties of membranes may further stabilize smaller do- 
mains in vivo. Recently, Yanagisawa et al. explored the 
consequences of a repulsive interaction between nearby 
buds [l2| and reported that such interactions can arrest 
the phase separation kinetics. The elastic perturbations 
induced by domains in the membrane, as described in this 
Letter, are obvious candidates for producing additional 
interactions between buds at any distance, further as- 
sisting in the creation of such a kinetic arrest. As Miiller 
et al. have shown for a flat membrane, two distortions 
on the same side _of an infinite flat membrane repel on 
all length scales 
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The experimental observation of 
multiple domains ordered in (quasi-)crystalline fashion 
in model membranes 14| strongly suggests a similar re- 
pulsive interaction in spherical vesicle systems. This is 
indeed evidenced by preliminary numerical exploration 
of this system using Surface Evolver [19]. It is of course 
straightforward to adapt the scheme outlined above to 
include long-range interactions between transmembrane 
proteins that impose a curvature on the membrane, e.g. 
scaffolding proteins [HI, lijj ]. Membrane mediated in- 
teractions act over length scales much larger than Van 
der Waals or electrostatic interactions and could provide 
an alternative or additional physical mechanism for pro- 
cesses like protein clustering and domain formation. Our 
results and methods allow not only to determine the pa- 
rameters relevant to processes like these, but also give a 
practical analytical handle on the shapes involved. This, 
in turn, will help decide between competing proposals for 
mechanisms involving membrane bending: protein inter- 
actions, endocytosis and the formation and stabilization 
of functional membrane domains. 
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